Elliptic flow from quark coalescence: mass ordering or quark number scaling? 



O 

o 

(N 
< 



(N 
> 

O 

00 

O 

o 



X 



Denes Molnar 

Department of Physics, Ohio State University, 174 West 18th Ave, Columbus, Ohio 43210, USA 

(Dated: February 8, 2008) 

We show that either mass ordering or quark number scaUng of anisotropic flow can result from 
quark coalescence, depending on the nature of phase space correlations at hadronization. Quark 
number scaling signals nonequilibrium dynamics because it can only appear when hydrodynamic 
correlations break down. Ifowever, the scaling does not hold for all nonthermal distributions, and 
is compatible with covariant transport theory only if remarkable cancellations occur at RHIC. 
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Introduction and conclusions. The goal of recent nu- 
clear collision experiments at the Relativistic Heavy Ion 
Collider (RHIC) is to create extremely hot and dense 
nuclear matter and determine its properties. An im- 
portant experimental probe in this quest is elliptic flow, 
V2 = (cos(20)), the second Fourier moment of the az- 
imuthal momentum distribution Its dependence on 
impact parameter b, transverse momentum, px and par- 
ticle species provides unique constraints on the equation 
of state (EOS) and collision dynamics. 

An amazing wealth of experimental V2{pt) data is 
available from Au + Au at ^snn = 130 and 200 GeV 
at RHIC for charged hadrons [1 E I3; 7r±, , and 
P + P |S@li7i"° 0; ^S' A + A 111; and + S+ and 
n- + n+ H. Below < 1 - 1-5 GeV, and for 5 < 8 
fm, these data can be reproduced fairl y w ell using ideal 
(Euler) hydrodynamics [iS [H E Hlli . In that the- 
ory, flavor dependence enters mainly through the mass 
of the particles, and the v^ipx) mass ordering pattern 
depends on the EOS [HlOEl- Remarkably, the data 
indicate a first-order deconfinement phase transition to 
the so called quark-gluon plasma (QGP) phase. 

In contrast, a spectacular quark number scaling [isllld 
H?! of elliptic flow V2 (pt) is predicted from hadronization 
via quark coalescence |3^9, 20, 21] . In the coalescence 
approach mesons (baryons) form from two (three) comov- 
ing quarks/antiquarks, leading to a unique meson-baryon 
differentiation in the flow pattern. 

Though mesons tend to be lighter than baryons, quark 
number scaling is in general different from mass ordering. 
Comparison between precise elliptic flow V2{pt) data for 
heavy mesons, such as the 0(1020) or K*(892), and light 
baryons (e.g., protons or lambdas) could distinguish be- 
tween the two behaviors. 

In view of the above difference, the realization that 
for thermal constituent distributions coalescence reduces 
to statistical hadronization generated signiflcant contro- 
versy. If, in that case, hadronization is the same both 
from coalescence and hydrodynamics, how could the re- 
sulting hadronic momentum anisotropics be different? 

In this paper we show that there is no contradic- 
tion at all. Either of the two scaling behaviors can re- 
sult from quark coalescence, depending on the type of 



phasespace correlations present at hadronization. Ther- 
mal constituent distributions imply unique coordinate- 
momentum correlations, for which mass ordering follows 
from coalescence. For quark number scaling to appear, 
departure from hydrodynamic behavior is essential. The 
growing experimental evidence 0,0,0 for quark number 
scaling in the intermediate transverse momentum region 
prp ^ 2 — 5 GeV at RHIC is a signal of the breakdown of 
hydrodynamics above pt > 2 GeV, which complements 
and corroborates other indications, such as the satura- 
tion of elliptic flow din m. 

Quark number scaling emerges only under certain con- 
ditions. Either only a subset of all possible local momen- 
tum anisotropics can be present, or important cancella- 
tions between contributions by the different anisotropics 
need to occur. Only the latter possibility is compatible 
with covariant transport theory, in which high-p^ parti- 
cles are emitted in a strongly preferred direction locally. 

Here we consider momentum anisotropics right after 
hadronization but before resonance decays. Secondary 
production affects signiflcantly the flnal anisotropics in 
hydrodynamics, but seems to influence little the quark 
number scaling pattern Q . At large ^ 5 GeV, contri- 
butions from jet fragmentation (ignored here) also need 
to be considered. 

In addition, this study considers coalescence on a 3D 
spacetime hypersurface. In coalescence from diffuse 4D 
freezeout distributions in spacetime, constituent space- 
time and space-momentum correlations influence the fi- 
nal hadron distributions in a more complex way [2ll | . 

Anisotropic flow. Momentum anisotropics can be 
characterized via the Fourier expansion of the freeze- 
out source distribution S{x,pT,y) = dN/d'^x d^pT dy in 
terms of the momentum azimuthal angle (j) 



S{x,pT,y) 

= So{x,pT,y) 



1 + 2 ^ Re {vn {x,pT,y)e~ 



-in(p\ 



(1) 



Here px = pt(cos 0, sin 0), (f> — is the impact param- 
eter direction, y = 0.51n[(i? +pz)/{E ~ Pz)] is the ra- 
pidity, and Vn{x,pT, y) € €' is the n-th order momentum 
anisotropy coefficient for particles emitted from an in- 



2 



finitesimal spacetime volume around x. The real and 
imaginary parts of u„ correspond to cos(n0) and sin(n(/i) 
terms. The final observable anisotropy u„ is the weighted 
average 



Vn{PT,y) 



J df-x Sn{x,PT,y) Vn{x,PT,y) 

J d'^x SQ{x,pT,y) 



(2) 



which is always real due to reflection symmetry across 
the collision plane (j) ~^ ~4' (for spherically symmetric 
nuclei). In addition, for a symmetric collision system, all 
odd Vn vanish at midrapidity y ~ due to — > (tt — 0) 
symmetry. 

In the following rapidity arguments will be often sup- 
pressed for brevity. 

Momentum anisotropies from ideal hydrodynamics. 
Two main assumptions in approaches based on ideal hy- 
drodynamics are local kinetic equilibrium and sudden 
freezeout on a 3D spacetime hypersurface. In such a case 

S{x,p) =p'^a^{x)5(x''cr^{x) - XoO-^(a;o)) f{x,p) (3) 

where cr^{x) is the normalized normal vector of the 
hypersurface. while xq is an arbitrary point on the 
hvpersurface |3 Ij . and 

where g is the degeneracy factor, and = — 1, 1, or 
respectively for bosons, fermions, or Boltzmann statis- 
tics. T, /i, and are the local temperature, chemical 
potential, and flow velocity. 

Momentum anisotropies have two sources, the hyper- 
surface (p ■ a) and the flow profile [p ■ u). For hyper- 
surfaces typically assumed in analytic studies, such as 
constant time or constant r = \fW^^z^ hypersurfaces, 
p • tr is independent of 4>. Therefore, for simplicity we fo- 
cus here on anisotropies due to the flow field only. Note, 
in general, 4> dependence can appear even for azimuthally 
symmetric hypersurfaces |25j . 

The anisotropy coefficients in general depend on the 
particle species, mainly because of the difference in mass. 
Statistics and the chemical potentials only play a role at 
very low mx = \/ m? + ~ (At freezeout at 

RHIC, T - 100 - 130 MeV, while ^ increases from ~ 80 
MeV for pions to ~ 350 MeV for heavy particles !l4i.) 

To demonstrate the mass dependence, consider not 
very low momenta, so that Boltzmann statistics is ap- 
plicable. The flow field can be characterized by a 
local transverse and longitudinal rapidity component, 
vt{x) = UT(cos$„,sin$„) and y{x), as m^(x) = 
(chy, WT, shy)/-\/l — v^. It is simple to show that 

V (x n^ v) - e™*"(") InishyTix)pT/T{x)) 
lo [shyT[x)pT/-L (x)) 



where shy^ = vt/ y/^ — v'x and {In} are the modified 
Bessel functions. Remarkably, the local anisotropy de- 
pends only on px, and therefore is the same for all par- 
ticles and rapidities. Note, at low pr, 



\vn\ « (shyTPT/2T)"/n! cx 



(6) 



as dictated by general analyticity arguments [2g, while 
at high pt, \vn\ ~ 1 — const x T/px — > 1. 

The coordinate-averaged flow coefficients on the 
other hand, in general decrease in magnitude with particle 
mass. This follows from the properties of the weight 5*0, 
which via ^ corresponds to the distribution 



/o = exp 



mTch(y - y{x))chyT{x) 
T{x) 



( pTshyxix) 



While \vn\ monotonically increases with the radial flow 
yr, the only mass dependent part in /o, the exponential, 
decreases with yx- The decrease is sharper for larger 
mass. Therefore for heavier particles, smaller Vn val- 
ues get preferred. Though this argument ignores the 
x-dependence of the cosine term from the exponential 
in ((Sl , which in principle could reverse the rise of Re w„ 
with yx and therefore affect the mass dependence pat- 
tern, the general trend prevails in practice. 

For more details of the derivations in the case of Ret;2, 
see H^, for example. 

Anisotropic flow from quark coalescence. In the quark 
coalescence approach, constituent quarks/antiquarks 
that are close in phasespace can combine to form hadrons. 
Assuming coalescence occurs on a 3D spacetime hy- 
persurface, the hadron binding energies are small, and 
coalescence is a relatively rare process 0, 0|, the 
invariant hadron momentum distribution can be ex- 
pressed 0, 0, 0, H^l in terms of the constituent phas- 
espace distributions and the hadron Wigner functions. 

Ignoring variations of phasespace distributions on 
length and momentum scales corresponding to a typical 
hadron {^-^ 1 fm and ~ 200 MeV), the phasespace distri- 
butions of mesons and baryons created via afi ^ M and 
a/37 B can be written as 



fB{x,p) 



(27r)3ga 
9a gf3 37 



fa{x,Pa.) fp{.X,pp)f^{x,p^) 



r / ^ (27r)6ffM , I ^ w / - N 
fM(X,p) = fa[X,Pa) fl3{x,Pi3) 

9a gp 



(7) 



Here g is the degeneracy of the particle (spin and color), 
while — P- Because constituents are comoving, 

p||pi and the hadron momentum is shared roughly in 
proportion to constituent mass [l^ . For hadrons com- 
posed of u, d, and s quarks, the sharing is approxi- 
mately equal because of the relatively small difference 
between niu^d ~ 0.3 GeV and rUg « 0.5 GeV. On the 
other hand, in hadrons that also contain a heavy quark, 
e.g., D mesons or the Ac, the heavy quark carries most 
of the momentum. 
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A direct implication of Q is 

a,/3,7 
i—a,f3 

where to leading order in Vn,i, the corrections Avn,M and 
Avn,B are 

oo n— 1 

Aw„ = ^ ■i;„+fe,iWfe + ^ v„-k,iVk,j + 0{{vi.i}^) 

i<j^k—l 

(9) 

with arguments identical to those in © dropped for 
brevity. Thus, if the corrections are small, w„ is additive, 
and in the absence of quark flavor dependence scales with 
quark number |l5lll6l |. 

Assuming, as for example in |l6| . that the spatial and 
momentum dependence of phasespace distributions fac- 
torize, the local Vn coefhcients are the same everywhere 
and therefore the scaling carries over to the observable 
spacetime averages Vn- In this case, the local Vn must of 
course be real, and at RHIC the odd ones must vanish at 
midrapidity. The neglect of the correction terms Av2 for 
elliptic flow is then justified because the data show [l^ 

Vq <^ V4 <^ V2 ^ 1- 

On the other hand, the scaling does not hold for ar- 
bitrary distributions. For example, for thermal (Boltz- 
mann) constituent phasespace distributions, coalescence 
(UI) gives thermal hadron distributions because, for the 
weak bound states assumed, mm ~ ma + mp, ms — 
ma + mf^s + m^ 1271] . Therefore, anisotropics are the same 
as from hydrodynamics, i.e., depend on the hadron mass. 
This also follows from a direct calculation. For mesons, 

oc 

the kernel is /o /o,a/o./3[l+ J2 (""".a^./s + c-c)], which 

71—1 

is identical to the hydrodynamic one with m = mM- The 
exponential part is the same because px = PT,a +PT,f3, 
momenta are shared such that my ~ TOT,a + 'mT,i3, and 
therefore also y ~ j/q — 2//3- The rest of the kernel and 
all meson flow anisotropics also agree, as can be shown 
using the exact flow addition formula 

Vn,M = 

oo n — 1 

i=a,(3 k=l ' k=l 

1 + E il'k,aVtB + CO.) 
fe=l 

(10) 

and the Bessel function addition theorem. The direct 
proof for baryons is analogous but more involved. 

Thus, for thermal constituent distributions, momen- 
tum anisotropics from quark coalescence depend on mass, 
and not quark number. There is an intuitive way to see 



why this must be so. In the frame where the flow ve- 
locity is zero, the momentum distribution is isotropic for 
all constituents. Therefore, in that frame all Vn vanish, 
and hence gives no hadron anisotropics. The only 
source of anisotropy then is the Lorentz boost back to 
the laboratory frame, which depends only on the particle 
mass, momentum, and the boost velocity. This is also 
true even if hydrodynamics breaks down, as long as dis- 
tributions are of the form f{x,p) — g{p ■ u{x), x). Quark 
number scaling can emerge only if, in any frame, at least 
one of the constituent distributions is anisotropic. 

In the thermal case, quark number scaling is violated 
because the nonlinear terms J^J are important. Consider 
for example V2, at low px so that is justified. To 
leading 0{p'^) order, both the linear and the vivi terms 
contribute. In fact the latter give as much as half the 
meson V2, and two-thirds for baryons. For higher-order 
flow coefficients, all Vn-k^k terms contribute at leading 
order in pT- 

Therefore, one class of freezeout distributions (besides 
uniform local anisotropy u„) for which the observable w„ 
scales with quark number is when 

\Vn\'>\Vn-k\\vk\ (11) 

in the spacetime region where most particles are emit- 
ted from. In this case, there is (approximate) scaling 
locally and, because the formula is linear, scaling is pre- 
served upon averaging. Note, a small amount of other 
anisotropics, for example vi Wj^^ besides a pure V2, 
can even help compensate the denominator in (|10|l , which 
otherwise tends to reduce the flow relative to the scaling 
expectation. 

The scaling may also hold even if the nonlinear terms 
are important. However, that requires a high degree 
of fortuitous cancellations at RHIC, as illustrated in 
Fig. ^ Results from the covariant parton transport 
model MFC [U IH are shown for Au + An at 
y/sNN ~ 200 GeV with 6 = 8 fm, from a calculation 
identical to the one in with a gluon-gluon cross sec- 
tion ffgg = 10 mb. In the left panel local constituent 
quark anisotropics up to fourth order are plotted as a 
function of transverse momentum, averaged over the first 
quadrant of the transverse plane (0 < tpx < 7r/2, where 
XT = XT (cos ipx, sin ifix))- Coefficients that vanish upon a 
full spacetime average due to symmetry, in particular vi 
and Imw2, are surprisingly large. Therefore, Hll() does not 
hold and, for example, the denominator in H10|l exceeds 
2 above pT,i > 1 GeV. 

The origin of these large anisotropics is shown in the 
right panel of Fig. ^ where we plot the (p distributions 
at midrapidity for 2 < < 3 GeV, averaged over four 
wedges in the transverse plane (px € [kir/S, {k + l)7r/8], 
with k — 0,1,2,3. Instead of small harmonic modula- 
tions over a uniform background, the distributions are 
strongly peaked because high-p^ particles can only es- 
cape from a surface layer of the reaction region. In this 
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FIG. 1: Results from MFC for Au + Au at ./a^ = 

200 GeV with 6 = 8 fm. Left: local constituent quark 
anisotropics at midrapidity as a function of pr, averaged over 
the first quadrant of the transverse plane. Right: (j) distribu- 
tions at midrapidity for 2 < pr < 3 GeV, averaged over four 
wedges in the transverse plane ipx G [kn/8, {k + l)7r/8], with 
k = 0...3. 



case, anisotropics from coalescence change their charac- 
ter completely. Consider, for simplicity, Gaussian con- 
stituent distributions around an azimuthal direction 0o, 
Si{x,pT,y) = a(x, pT, 2/) c xp[-(0 - <Po)^/{2af)]. For 
not too wide a < 2/y/n + 1, the constituent and hadron 
anisotropics 

(12) 

scale as {const)~'^ . This is very different from both 
and For example, if all constituents have the same 

width cr, we have 

k«,s(3pT)| ^ k„,g(pT)^/^ b„,M(2pT)| ^ \Vn^gipT)\^/^ , 

which gives for ^0.1 smaller baryon/meson flow ra- 
tios, and for \vn^q\ > 0.25 also smaller hadron flows, than 
Vn,s(3pT) = 3w„,g(pT), Vn,M{2pT) = 2w„^,(pt) from lin- 
ear scaling (0). 

The observable averages are, of course, deter- 

mined by the interplay between variations in the local 
emission angle, the width, and constituent density. It 
is quite remarkable that the end result in a dynamical 
coalescence approach [2^ is only a modest ~ 20% and 
^ 30% reduction of pion and proton elliptic flow at RHIC 
relative to quark number scaling. 

The above results underscore the importance of local 
momentum anisotropics in hadronization via quark co- 
alescence. Flow addition formulas apply locally, even 
though only spacetime averages of the final anisotropics 
can be observed. Odd-order and the sine terms 

play a major role, even at midrapidity. These terms have 
been ignored in analytic studies so far 0, js^l but are 
naturally incorporated in a dynamical approach 

Acknowledgments. Helpful discussions with E. Shuryak 
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